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Abstract. Tube formulas (by which we mean an expHcit formula for the vol- 
ume of an £-neighbourhood of a subset of a suitable metric space) have been 
used in many situations to study properties of the subset. For smooth subman- 
ifolds of Euclidean space, this includes Weyl's celebrated results on spectral 
asymptotics, and the subsequent relation between curvature and spectrum. 
Additionally, a tube formula contains information about the dimension and 
measurability of rough sets. In convex geometry, the tube formula of a con- 
vex subset of Euclidean space allows for the definition of certain curvature 
measures. These measures describe the curvature of sets which are not too 
irregular to support derivatives. In this survey paper, we describe some recent 
advances in the development of tube formulas for self-similar fractals, and their 
applications and connections to the other topics mentioned here. 



1. Introduction 

This survey article describes advances on the computation of tube formulas for 
fractal subsets of M and of M"*, and relations to classical results. In particular, we 
show how the theory of complex dimensions can be used to calculate explicit tube 
formulas for a large class of self-similar fractals. We also discuss generalizations 
and connections to forthcoming work. 

We begin by discussing results for fractal strings, that is, fractal subsets of R. 
Also, we give the generalization to tube formulas of measures, and show how these 
results may be applied to the investigation of dimension and measurability of rough 
sets. In particular, we use the central notion of "complex dimensions" , a concept 
which extends the real-valued notion of Minkowski dimension (and Hausdorff di- 
mension, in some cases). §3 describes the explicit computation of a tube formula 
for the von Koch snowflake curve, a classical subset of with fractal boundary. 
This computation is mainly a "brute-force" calculation using elementary geometry 
coupled with some subtle exploitation of distributional methods. §4 contains the 
construction of a self-similar tiling via an iterated function system, and shows how 
the tiling enables one to extend the results for fractal strings to self-similar subsets 
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of R . §5 gives a brief description of some classical results from convex geome- 
try and shows how these relate to the results for self-similar tilings. Finally, we 
conclude with a pair of contrasting examples that illustrate the key ideas discussed. 

1.1. Basic concepts. We now present some of the ideas that are used through- 
out, and introduce basic notation. 

Definition 1.1. Given e > 0, the inner e -neighbourhood of a set ^ C M'*, 
d > 1, is 

As :={xe A ; dist{x, dA) < e}, (1.1) 

where dA is the boundary of A. For a given A, we are primarily interested in 
finding a tube formula for A, that is, an explicit expression for the d-dimensional 
Lebesgue measure of A^, denoted 

Va{s) -.^voUiA,). 

Definition 1.2. A self-similar system is a family (with J > 2) of 

contraction similitudes 

^j{x) := tjAjX + aj, j = 1, . . . , J. 

For j = 1,...,J, we have < < 1, aj G K'^, and Aj G 0{d), the orthogonal 
group of rigid rotations in d-dimensional Euclidean space M"*. The number tj is the 
scaling ratio of $j. For convenience, assume that 

l>ri>r2>--->rj>0. (1.2) 

When d = 1, one has only Aj = ± 1 . 

It is well known that there is a unique nonempty compact subset F C M.'^ 
satisfying the fixed-point equation 

J 

F = $(F):= y $,(F). (1.3) 

This (self-similar) set F is called the attractor of Given a word w = wi . . .Wk S 
Wk := {1, 2, . . . , J}*^, wc denote the composition of several similarity mappings by 

2. Fractal strings 

The essential strategy of fractal strings is to study fractal subsets of K by 
studying their complements. 

Definition 2.1. A fractal string is any bounded open subset L CM, that is, 
a countable sequence of open intervals 

L :={i„}^=i. (2.1) 

Although it is not part of the definition, the idea is that the boundary dL 
is a set one wants to study. When analyzing L in terms of characteristics which 
are invariant with respect to rigid motions, the position of an interval L„ within 
R is immaterial, and all pertinent data is stored in the sequence of lengths £n of 
Ln- Indeed, this is the formulation in which the concept is used in [La-vF2], the 
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primary reference for fractal strings. For agreement with the sequel, however, we 
use the inradius. 

Definition 2.2. The inradius p of a set A is 

p = p[A) := sup{e > ; 3a; with B{x, e) C A}. (2.2) 

It is clear that if ^ is a bounded set, AC A^ for sufficiently large e. Alterna- 
tively, it is apparent that for a fixed e > 0, any sufficiently small set will be entirely 
contained within its £-neighbourhood. The notion of inradius allows us to see when 
this phenomenon occurs; its relevance should be clear from Definition 1.1. 

Obviously, the inradius of an interval is p„ = = 5Voli(i„), i.e., half the 
length of the interval L„ C M. We then divide by the first (largest) length to 
normalize the inradii. The result is a sequence of scaling ratios (or scales) {r„}, 
and ii g is the length of the largest interval, then = 2grn. Use of the terminology 
"scaling ratio" here corresponds to the implicit idea that each interval Ln £ L 
is congruent to a copy of the largest interval which has been scaled by r„. For 
a self-similar string, this mapping is explicitly given by <!>„,; see (2.4). Thus, by 
the remarks following Definition 2.1, all pertinent data is stored in the sequence of 
scaling ratios 

C :={r„}- 1. (2.3) 

We may take the sequence in nonincreasing order so that ri > r2 > • • • > 0. To 
avoid trivialities, assume that r„ > for all n e N. Of course, since L is bounded 
we have X^^i < 

To define a self-similar string, consider the self-similar system with similarity 
mappings ^j{x) = tjAjx + aj. The set of scaling ratios of a self-similar string will 
consists of the collection of all products of scaling ratios of the maps. In particular, 
if w is a finite word on {1, 2, . . . , J}, then the string will include 

the scaling ratio of = ^w^. o . . . o ^yj^ o ^yj^ . Thus, a self-similar string contains 
every number that arises as the scaling ratio of a composition of the similarity 
transformations . (The first one is always ri = 1 , corresponding to a composition 
of similarities.) The motivation for this definition is that it corresponds to the 
set obtained by taking the set-theoretic difference of a self-similar subset of R from 
the smallest closed interval containing it. Alternatively, the set may be constructed 
by selecting any bounded open interval ICR and examining the lengths of the 
intervals {^^{I)}, where w runs over all finite words. An example of a self-similar 
string is given by the Cantor string in Example 2.8 just below. 

A generalized fractal string is a locally finite Borel measure on (e, oo), where 
e > 0, and is denoted r] = rjc- Such a string may not have a geometric realiza- 
tion. The motivation for this generalization lies in the flexibility of working in the 
measure-theoretic framework, and in certain specific applications. In this context, 
an ordinary fractal string as given in the previous definitions (self-similar or not) 
corresponds to a sum of Dirac masses Sx, each located at a reciprocal of one of the 
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numbers r„. That is, an ordinary fractal string may be written 

oo 

% = ^V„- (2-5) 

Definition 2.3. The scaling zeta function (s of a fractal string is the Mellin 
transform of the measure r]c- 

POO 

Us) := / x-'drjcix), see. (2.6) 
Jo 

In [La-vF2], this is called the geometric zeta function of C We have chosen the 
current terminology to agree with the latter sections of this paper and instead say 
that the geometric zeta function of (c is given by 

Cc{e,s):=Us)^^y (2.7) 

The function (c factors into and a term which contains geometric data about 
open intervals, although this will not become clear until the discussion of the higher- 
dimensional tilings in §4. In the case when C is self-similar, the function : f2 — > C 
may be meromorphically continued to all of C. Otherwise, it is well defined in some 
half-plane {Re s > a} and may be continued analytically. To demarcate the domain 
of (c formally, we introduce the following definition. 

Definition 2.4. Let / : M — > M be a bounded Lipschitz continuous function. 
Then the screen is S = {f{t)+it t e M}, the graph of / with the axes interchanged. 
The region to the right of the screen is the window W := {z G C \ Re2;> 
/(Imz)}. We choose / so that S does not pass through any poles of (c, and (c has 
a meromorphic extension to some neighbourhood of W. 

Definition 2.5. The (complex) scaling dimensions of C arc poles of <^£. The 
poles which lie in the window are called the visible scaling dimensions: 

Vc{W)={ujeW I lim\Ccis)\=oo}. (2.8) 

The screen and window are useful for many purposes. In addition to demarcat- 
ing the domain of , they allow one to make precise statements about the growth 

of (cy and they allow for certain approximation arguments that make the proof of 
the next theorem possible. In particular, one has the following definition. 

Definition 2.6. One says that (c (or just £) is languid if it satisfies certain 
mild growth conditions on S, and along a sequence of horizontal lines in W. For 
the precise nature of these conditions, see [La-vF2, Def. 5.2]. 

Theorem 2.7 (Tube formula for fractal strings [La-vF2], Thm. 8.1). Let rjc 
be a fractal string with geometric zeta function C,c and assume that (^c is languid. 
Then we have a tube formula 

Vc{e) = ^ res {Ccie, s); s = w) + {2eC^(0)} + TZ{e). (2.9) 

Here the term in braces is only included ifO&W\ ^^{W), and the error term is 

n{e) = ^ [ Cc{£,s)ds = 0{e^-^''P^), ase^O+. (2.10) 
27ri Js 

When the string is also self-similar, one may take W = C and TZ{s) = 0. 
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If we denote the poles of the scaling zeta function separately by , then the 
result above may be rewritten 

Vc{s)= Yl res(C£(e,s);s = a;)+7e(e). (2.11) 

Formula (2.9) was originally obtained as a distribution acting on smooth functions 
with compact support in (0, oo). However, it has since been obtained in a pointwise 
fashion in [La-vF2, Thm. 8.7] under only slightly more restrictive conditions (which 
are always satisfied in the case of self-similar strings), so this technicality need not 
be emphasized. 

Example 2.8. The complement of the usual Cantor set in the unit interval 
[0,1] consists of open intervals with lengths This self-similar 

fractal string has a largest interval of length i and scaling ratios 3"*^ appearing 
with multiplicity 2*^. Consequently, the Cantor string CS may be written 

C30 

%5 = 5Z2%., (2.12) 

fe=0 

and its scaling zeta function is 

Us) = l ^-'drics{x)=^^i^-j = ^_^.3_, , (2.13) 

while its tube formula is 

1/1 1 \ /£\i~^~'"P 

with g = I is the inradius of the largest interval, D = log^ 2 is the Minkowski 
dimension, and p = 27r/log3 is a constant called the oscillatory period. 

Definition 2.9. The Minkowski dimension of A C R'^ is 

D = dimM A := mi{t > ; ^^(e) = 0{e'^-*), as £ ^ 0+}, (2.15) 

and is also frequently called the box dimension. For a string, dim£ := dim9L. 

The complex dimensions can be thought of as a generalization of Minkowski 
dimension because of the following result of [La2], which also appears as [La- 
vF2, Thm. 1.10]. 

Theorem 2.10. For a string £, the scaling zeta function Cs converges on the 
half-plane {s e C ] Res > a} if and only if a > dim£. 

Remark 2.11. In formula (2.13), one discovers the pleasant surprise that has 
a nice closed form. In fact, it is shown in [La-vF2, Thm. 2.4] that all self-similar 
strings have a scaling zeta function of the form 

Us) = -^^. (2.16) 

This remark remains true for the self-similar tilings discussed in §4. The number- 
theoretic and measure-theoretic implications of this result are extensive; they are 
studied at length in [La-vF2, Ch. 2-3]. 
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Figure 1. The Koch curve K and Koch snowflake domain O. 




Figure 2. Four early stages in the geometric construction of K. 



3. The von Koch curve 

The Koch curve K is the attractor for the self-similar system 

$i(z):=pz and $2(2) (1 - - 1) + 1, (3.1) 

where p = ^ -I- and i = Consequently, K is & self-similar fractal with 

Minkowski dimension D log3 4; see Figure 1. In this section, we describe how 
the tube formula for the Koch snowflake was calculated directly. 

hi [LaPel], the authors compute a tube formula for the Koch snowflake Q. 
by considering a sequence of curves Kn —f K, where convergence is with respect 
to Hausdorff metric. Some representative terms of this sequence are illustrated in 
Figure 2. The tube formula is obtained as a limit of tube formulas Vk„ obtained by 
approximation for each K^, as illustrated for K2 in Figure 3. Three copies of this 
figure are fitted together to form the inner neighbourhood of the snowflake, as is 
indicated by the dashed lines at either end. This region is only an approximation of 
course, as one side of each rectangle should be replaced with a fractal curve. Fig. 4 
shows how this error is incurred and how it inherits a Cantoresque structure from 
the triadic character of the Koch curve. Let us refer to the dark region in Fig. 4 
as an error block and each connected component Ak of it as a trianglet. Without 
taking the error blocks into account, the e-neighbourhood of the Koch curve has 
approximate area 

Vk(£) = £2-C4-M (3^9{^> + ^3{^> + 1 (1 _ Vs)) - 4 (f + 2V3) . (3.2) 

Here, x := — \og^{e^/i) , and we use x ~ [x] + {x} to denote the decomposition of 
X into its integer and fractional parts. That is, [x] is an integer and < {x} < 1. 
Furthermore, x is related to n (the index of Kn) by n = n(e) = [log3 = [x\. 
Thus, the level n of the approximation is determined by e, with n — > 00 as £ ^ 0. 



TUBE FORMULAS FOR SELF-SIMILAR FRACTALS 



7 




Figure 3. An approximation to the inner e-neighbourhood of the 
Koch curve, based on the curve K2 from Figure 2. 




Figure 4. An error block for Kn- The central third of the block 
contains one large isosceles triangle, two wedges, and the trianglet 
Ai. Figure 5 contains 4 complete copies of this figure, and 12 
partial copies of it. 

The tube formula for il will be obtained by summing the areas Ak of the 
trianglets, and multiplying by the number of error blocks occurring for a given 
approximation, i.e., for a particular value of e. Unfortunately, this number is not 
easy to express, due to the existence of partial error blocks. See Figure 5 for a 
visual explanation of what is meant by partial error blocks and complete error 
blocks. The number of complete error blocks can be readily counted with a simple 
formula, but the portion of a partial error block that exists for a given value of 
e is rather ornery, and so we denote this quantity by h{e). We do not know h[e) 
explicitly, but we do know by the self-similarity of K that it has multiplicative 
period 3; i.e., h{e) = It is a theorem of [LlWi] that for a self-similar set, the 

e-neighbourhood must be either: (i) a finite union of convex sets for every value of 
s, or (ii) not a finite union of convex sets for any value of e. We conjecture that for 
every self-similar set of the latter type, the tube formula involves a multiplicatively 
periodic function analogous to our h{s). Once the error has been taken into account, 
one obtains [LaPel, Thm. 5.1] as follows. 
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partial 




Figure 5. Error block formation. The ends are counted as partial 
because three copies of this illustration will be added together to 
make the entire snowflake. 




J !_ J_ _l !_ ± 



9V5 SVS Vi 9\/3 3V5 v5 

Figure 6. A comparison between the graph of the Cantor- like 
function h and the graph of an approximation h{e) = ^ - {—[x]—x), 
where /x is a constant and x = — log3(£\/3)- 

Theorem 3.1. The area of the inner e -neighbourhood of the Koch snowflake is 
given by the following tube formula: 

Vis) = J2 V'r.e'-^-^P + V'ne'""^ (3-3) 

for suitable constants tpn , ■0n which depend only on n and are expressed in terms of 
the Fourier coefficients g^ of h{e). 

To see the full form of (3.3), please see Remark 6.4 at the end of this paper. 
The result (3.3) was obtained at a time when the theory of complex dimensions 
was entirely restricted to fractal strings as outlined in §2 (except for the conjectures 
expressed in [La-vF2, Ch. 12]). Part of the motivation for proving Theorem 3.1 
was to get an idea of what the complex dimensions of Q, might look like. Reasoning 
by analogy, one would deduce from (3.3) that the complex dimensions of 11 are 
obtained from the exponents appearing in (3.3); in particular, that each exponent 
is of the form 2 — w, where w is a complex dimension of O. This led to the prediction 
that the complex dimensions of the Koch snowflake are 

'Ddn = {D + inp ; n G Z} U {inp ; n e Z}. (3.4) 

The reader will find in Example 6.2 and Remark 6.4 that this is not far wrong. 
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4. Fractal sprays and self-similar tilings 

A fractal spray is the higher-dimensional counterpart of a fractal string. See 
[La-vF2, § 1.4] for a discussion; this idea also appears earlier in [LaPo2], [La2], 
and [La3]. 

Definition 4.1. Let G C K'' be a nonempty bounded open set, which we 
will call the generator, and let C = {r„} be a fractal string. Then a fractal spray 
is a bounded open subset of R'^ which is the disjoint union of open sets i?„ for 
n = 1,2, . . . , where each i?„ is congruent to r„G, the homothetic of G by r„. 

Thus, any fractal string can be thought of as a fractal spray on the basic shape 
G = (0, 1), the unit interval. Every self-similar system is naturally associated to a 
certain fractal spray called the self-similar tiling. 

Definition 4.2. The self- similar tiling T corresponding to a self-similar system 
$ is a fractal spray where the fractal string and generators are defined as follows. 
The string is the collection of all finite products of the scaling ratios {ti, . . . , tj} of 
the self-similar system $. Let C = [F] be the convex hull of F, and let int C be its 
interior. Then the generators are the connected open sets Gg in the disjoint union 

intG\$(G) =GiU---UGq. (4.1) 

When there is more than one generator, it is more accurate to think of the 
tiling as a union of fractal sprays, one for each generator. Indeed, for purposes 
of computing the tube formula, it is easiest to deal with each part separately and 
then obtain the final result by adding the contributions from each generator. For 
this reason, the tube formulas below are all phrased for a spray or tiling with one 
generator. In fact, we cannot currently exclude the possibility that there may exist 
examples for which Q = co. However, we have been unable to construct such an 
example, and for the time being we assume Q < oo. 

The term "self-similar tiling" is used here in a sense quite different from the one 
often encountered in the literature. In particular, the tiles themselves are neither 
self-similar nor are they all of the same size; in fact, the tiles are typically simple 
polyhedra. Moreover, the region being tiled is the complement of the self-similar 
set F within its convex hull, rather than all of M.'^; see Figure 7. In fact, it is shown 
in [Pe2] that when $ satisfies the tileset condition, the collection {^w{Gq)}w,g 
forms an open tiling of C\F. We now define these terms. 

Definition 4.3. An open tiling of A gW^ is a, collection of nonempty connected 
open sets {An}^^i such that 

] [oo 

{i)A = [\ An, and (ii) A„ n = forn ^ m. (4.2) 

Figure 7 shows a tiling by open sets which are equilateral triangles. 

Definition 4.4. The system $ satisfies the tileset condition iff T ^ ^(C) and 

int$j(G)nint$£(G) =0, j £. (4.3) 

The tileset condition is a separation condition which is similar to the "open set 
condition", but stronger. The open set condition is satisfied when there exists a 
"feasible open set" U which has the property that ^j{U) C U, and ^j{U) n ^e{U) 
for all j, £. To see that the tileset condition implies this, let U be the the interior 
of G = [F]. For a counterexample to the converse, see [Pe2, Example 3.8]. 
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For self-similar tilings, and more generally for fractal sprays, the scaling zeta 
function is defined just as it was for fractal strings. However, the geometric zeta 
function becomes much more complicated, due to the multifarious possibilities for 
the geometry of the generator. The resulting technicalities can be ignored in many 
cases, however, for example when the generator is sufficiently simple. This mo- 
tivates the condition "Steiner-like" in the definition just below. A bounded open 
set in M.'^ is said to be Steiner-like if its inner parallel volume Vg{s) admits a 
"polynomial-like" expansion in e of degree at most d. More precisely, we have the 
following definition (see §5 for an explanation of this choice of terminology). 

Definition 4.5. A bounded open set G C M^' is Steiner-like iff for < e < p{G) 
its inner tube formula may be written 

d 

VG{e)=J2^k{G,e)e''-\ (4.4) 

where each coefficient function Kk{G,-) is assumed to be a bounded and locally 
integrable function of e for which 

lim Kfc(G,e) (4.5) 

exists, and is both positive and finite. 

Definition 4.6. In the special case when G is a Steiner-like set whose terms 
are constant, we say the set is diphase. This terminology refers to the fact that its 
tube formula is written piecewise with only two cases: e < p{G) and e > p{G). If G 
is Steiner-like and the functions Kfe are piecewise constant on the interval [0,p(G)], 
then we say G is pluriphase. Thus, diphase is a special case of pluriphase. 

Conjecture 1. The class of pluriphase sets includes convex sets and polyhe- 

dra. 

This result may exist in the literature, but we have been unable to find it. 
It is simple to show that the class of diphase sets includes spheres and regular 
polyhedra. In particular, the examples of §6 are both diphase. There are examples 
of other types of diphase sets, but they resist easy description. Current research is 
attempting to characterize diphase and pluriphase sets in terms of other geometric 
properties. 
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In the case when G is pluriphase, the geometric zeta function is relatively simple 
to write down. To avoid obscuring the exposition, we give only the diphase case 
here (although in view of Definition 4.8 it is worth noting that the factors also 
appear in the more general case). 

Definition 4.7. The geometric zeta function of a fractal spray with a diphase 
generator is 

Cr (e, s) ■.=s''-%r,{s) (4-6) 

Definition 4.8. The set of visible complex dimensions of a fractal spray is 

Vr{W) := Vr,{W) U {0, 1, . . . , - 1}. (4.7) 

Thus, VriW) consists of the visible scaling dimensions and the "integer dimen- 
sions" of G, and contains all the poles of Qr (viewed as a function of s). 

Theorem 4.9 (Tube formula for fractal sprays). Let rj be a fractal spray on the 
Steiner-like generator G, with generating inradius g = p{G) > 0. If (t is languid, 
then we have a tube formula 

Vr{e)= res(Cr(e,s);a;)+7^(e), (4.8) 

where the sum ranges over the set VriW). Here, the error term 'R-{£) is 

TZ{e) = ^ [ Crie, s) ds, = 0(£'^-^"p^), as e ^ 0+. (4.9) 
27ri Jg 

Theorem 4.9 was first obtained distributionally in [Pel] and appears in im- 
proved form in [LaPe2]. Since then, a pointwise proof has been obtained in 
[LaPeWil]. It is important to note that Cr is languid for all self-similar tilings. 
As we have not been able to construct a self-similar tiling which satisfies the tileset 
condition but fails to be Steiner-like, Theorem 4.9 applies to all known examples. 
A fortiori, it is possible to show that all self-similar tilings automatically satisfy a 
more stringent condition (called strongly languid in [La-vF2, Def 5.2]) that enables 
one to take W = C and TZ{s) = 0; see [LaPe2, Thm. 8.4]. Thusly one obtains the 
following special case of Theorem 4.9. 

Corollary 4.10. Let T be a self-similar tiling with a diphase generator G. If 
Ct{s) has only simple poles, then 

d d-l 

^^(^)= E E^^«(C.(s);a;)£^-^«, + ^«,a(^)£'-^ (4.10) 

w£l>, 1=0 i=0 

By comparing (4.10) with (2.9), it is easy to see how the tube formula for fractal 

sprays extends the results for fractal strings. Additionally, it extends classical 
results for convex sets, as outlined in the next section; see especially (5.6)-(5.7) 
and the surrounding discussion. 
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5. Convex geometry and the curvature measures 

In order to explain the connections with convex geometry, we give a brief en- 
capsulation of Steiner's classical result. Here, we denote the Minkowski sum of two 
sets in by 

A + B ^ {x e ■ X = a + b for a e A,b e B}. 

Theorem 5.1. If B"^ is the d- dimensional unit hall and A^W^ is convex, then 
the d- dimensional volume of A + eB'^ is given by 

d 

YoldiA + sB") = ^ ^^i{A) yold.i{B''-')e''-\ (5.1) 

i=0 

where fXi is the renormalized i-dimensional intrinsic volume and Yo\j{A} is the j- 
dimensional Lebesgue measure. 

Up to some normalizing constant, the z-dimensional intrinsic volume is the same 
thing as the i^^ total curvature or (d — i)**^ Qucrmassintegral. This valuation /Zj can 
be defined via integral geometry as the average measure of orthogonal projections 
to [d — i)-dimensional subspaces; see Chap. 7 of [KlRo]. For now, we note that 
(up to a constant), there is a correspondence 

jjbQ <^ Euler characteristic, Md-i ~ surface area, 

Hi mean width, fi^ volume, 

see [Schn2, §4.2] for more. We have chosen the term "Steiner-like" for Def. 4.5 
because the intrinsic volumes satisfy the following properties: 

(i) each /Uj is homogeneous of degree i, so that for a; > 0, 

{xA) = iJLi{A) x\ and (5.2) 

(ii) each \Xi{^A) is rigid motion invariant, so that 

Mi (T(A)) = ^x^{A), (5.3) 

for any (afRne) isometry T of M''. 

Note that (5.1) gives the volume of the set of points which are within e of A, 
including the points of A. If we denote the exterior e -neighbourhood of A by 

* := {A + eB'^) \A = {x \ d{x, A)<e,x^ A}, 
then it is immediately clear that omitting the d*"^ term gives 

d-1 

vold(^r*) = 5^Ci(A)e<*-' (5.4) 

i=Q 

with Ci{A) = Mj(A) vold.i(i?'^~*). The intrinsic volumes fXi can be localized and 
understood as the curvature measures introduced in [Fed] and described in Ch. 4 
of [Schn2]. In this case, for a Borel set (3 C R*^, one has 

d-1 

yoU{x G Af * ; p{x,A) e (3} = (A, /?)£''"', (5.5) 

where p{x, A) is the metric projection of x to A, that is, the closest point of A to 
x. In fact, the curvature measures are obtained axiomatically in [Schn2] as the 
coefficients of the tube formula, and it is this approach that we hope to emulate 
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in our current work. In other words, we believe that may also be understood 
as a (total) curvature, in a suitable sense, and we expect that m can be localized 
as a curvature measure. A more rigorous formulation of these ideas is currently 
underway in [LaPeWiS]. Caveat: the description of Ki given in the conditions of 
Def. 4.5 is intended to emphasize the resemblance between Ki and Cj. However, Ki 
may be signed (even when G is convex and i = d — l,d) and is more complicated 
in general. In contrast, the curvature measures Cj are always positive for convex 
bodies. 

The primary reason we have worked with the inner £-neighbourhood instead 
of the exterior is that it is more intrinsic to the set; it makes the computation 
independent of the embedding of T into W^. At least, this should be the case, 
provided the 'curvature' terms Ki of Def. 4.6 are also intrinsic. As a practical 
bonus, working with the inner e-neighbourhood allows us to avoid potential issues 
with the intersections of the e-neighbourhoods of different components. 

In [Fed], Federer unified the tube formulas of Steiner (for convex bodies, as 
described in Ch. 4 of [Schn2]) and of Weyl (for smooth submanifolds, as described 
in [BeGo], [Gr] and [We]) and extended these results to sets of positive reach? It is 
worth noting that Weyl's tube formula for smooth submanifolds of is expressed 
as a polynomial in e with coefficients defined in terms of curvatures (in a classical 
sense) that are intrinsic to the submanifold [We]. See §6.6-6.9 of [BeGo] and the 
book [Gr]. Federer's tube formula has since been extended in various directions 
by a number of researchers in integral geometry and geometric measure theory, 
including [Schnl,Schn2], [Zal,Za2], [Ful,Fu2], [Sta], and most recently (and 
most generally) in [HuLaWe]. The books [Gr] and [Schn2] contain extensive 
endnotes with further information and many other references. 

To emphasize the present analogy, consider that Steiner's formula (5.4) may be 
rewritten 

vold(Ar*)= ""i^"''- (5-6) 

ie{0,l,...,d-l} 

and it is clear that Corollary 4.10 may be rewritten 

Vr{e)= c-^''""' (5.7) 

a'eI>,U{0,l,...,d-l} 

where for each fixed w G P^, 

Ca, := res(C5(s);w) V -.Ki. (5.8) 

Note that when a; = i G {0, 1, . . . , d — 1}, one has c^j = Cj = ^^(z)^^. The obvious 
similarities between the tube formulas (5.6) and (5.7) is striking. Our tube formula 
is a fractal power series in e, rather than just a polynomial in e (as in Steiner's 
formula). Moreover, our series is summed not just over the 'integral dimensions' 
{0, 1, . . . , d — 1}, but also over the countable set of scaling complex dimensions. 
The coefficients c^j of the tube formula are expressed in terms of the 'curvatures' 
and the inradii of the generators of the tiling. 

set A has positive reach iff there is some S > such that any point x within 5 of A 
has a unique metric projection to A, i.e., that there is a unique point A niiiiiinizing dist{x, A). 
Equivalently, every point q on the boundary of A lies on a sphere of radius <5 which intersects dA 
only at q. 
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A 




A 




Figure 8. The Sierpinski gasket tiling. The first stages of the 
construction of the gasket are labeled C„. The tiles {^^{G)} are 
labeled r„ where — n. 

Remark 5.2. The two formulas (5.6) and (5.7) initially appear to be measuring 
very different things, but this may be misleading. It can happen that the exterior 
£-neighbourhood of the fractal itself is, in fact, equal to the union of the inner 
e-neighbourhood of the tiling and the exterior e-neighbourhood of its convex hull. 
In this case, (5.6) and (5.7) are closely related. In fact, the tube formula for the 
exterior e-neighbourhood of the fractal is precisely the sum of the tube formula for 
the tiling and Steiner's tube formula for its convex hull. This occurs, for example, 
with the Sierpinski gasket; see Remark 6.5 and Figure 10. 

Remark 5.3 (Comparison of Vq- with the Steincr formula). In the trivial sit- 
uation when the spray consists only of finitely many scaled copies of a diphase 
generator (so the scaling measure is supported on a finite set), the scaling zeta 
function will have no poles in C, so that 7)^ = 0- Therefore, the tube formula 
becomes a sum over only the numbers i = 0, 1, . . . , — 1,^ for which the residues 
simphfy greatly. In fact, in this case C?;(«) — Pi + ' ' ' + Pjvj so each residue from 
(4.10) is a finite sum 



C^(i)K,(e) = p\n,e'^-' + ■■■ + PNK^e'^~' = K,{r^,G)e'^-' + ■■■ + K,{r,,^G)e'^-' 



where N is the number of scaled copies of the generator G, and is the corre- 
sponding scaling factor. Summing over i, we obtain a tube formula for the scaled 
basic shape r^^G, for each n— 1, . . . ,N. The pluriphase case is analogous. 



Example 6.1 (The Sierpinski gasket tiling). The Sierpinski gasket tiling SQ 
(see Fig. 8) is constructed via the system 



which has one common scaling ratio r = 1/2, with J = 3 and one generator G: 



6. Two illustrative examples 



$i(z):=iz, <i>2{z) := ^z + ^, <S>3{z) := + 



l+iV3 



4 




■^It turns out during the computation that the d^^ summand always vanishes, so the sum 
extends only up to d — 1. 
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Figure 9. A measurable Koch tiling obtained by picking a non- 
lattice parameter ^ in (3.1). 

and the scaling complex dimensions are 

V^^{D + inp \ne1} for Z) = log2 3, p = (6.2) 
The tube formula for SQ is readily computed to be 

V (A = \^ ( l_ , _^ \ (e\ , t!±^2 _ 

^Sg\t.j l&\og2 \ D+inp ' Z5-l+inp D-2+inp ) \g ) 2 ^ 

Example 6.2 (The Koch tiling). The standard Koch tihng /C (see Figure 7) is 
constructed via the self-similar system given in (3.1). The attractor of {$i,$2} is 
the classical von Koch curve as in Figure 1. This system has one scaling ratio 
r = 1^1 = l/-\/3, with J — 2 and one generator G: an equilateral triangle of side 
length i and generating inradius g — -j^- This tiling has inradii pk ~ gt^ with 
multiplicity 2'^, so the scaling zeta function is 

QJs) = ^- ^, (6.3) 

and the scaling complex dimensions are 

= {i? + mp;neZ} for D = logg 4, p = (6.4) 
By inspection, a tile with inradius 1 /x will have tube formula 

/33/2(-£2 + 2ex), £<l/x, 
^^^^'^^ = \33/^.^ e>l/x. 

For fixed x, (6.5) is clearly continuous at e = 0+. Thus we have kq = — 3^^^, ki — 
2 • 3^/^, and K2 = — 3^/^. Now applying (4.10), the tube formula for the Koch tiling 
JC is 

^'^^^^ ^ bg3 ^ D+\np + _D-l+i«p ^ Ll-2+inp) (f ) 



+ 33/V+ ^_^.1_,,, £, (6.6) 



where D = logg 4, 5 = ^ and p as before. 

Remark 6.3 (Nonlattice Koch tilings). By replacing ^ = 5 + in (3.1) with 
any other complex number satisfying + < one obtains a family of 

examples of nonlattice self-similar tilings, as illustrated in Figure 9. Computation 
of the tube formula does not differ significantly from the lattice case. Further dis- 
cussion of nonlattice Koch tilings may be found in [Pe2] . As discussed in [LaPe2] , 
for example, this furnishes a family of tilings, almost all of which are Minkowski 
measurable. 



16 



M. L. LAPIDUS AND E. P. J. PEARSE 



Remark 6.4. In §3 wc discussed how to obtain a tube formula was obtained for 
the £:-neighbourhood of the Koch curve itself (rather than of the tiling associated 
with it) and this led to the prediction that the complex dimensions of the curve are 

VtCi, = {D + inp : n e Z} U {0 + inp ; n G Z}, 

where D = \og^4, and p = j^g. The line of poles above D was expected^, and 
agrees precisely with the results of §4. The meaning of the line of poles above is 
still unclear. Wc invite the reader to compare (6.6) with the formula for the Koch 
curve. When formula (3.3) for the area of the inner e-neighbourhood of the Koch 
snowflake is stated in full detail, it appears as follows. 

V{e) = Gi(e)£'-^ + G'2(£)£', (6.7) 

where G\ and G2 are periodic functions of multiplicative period 3, given by 

Gi(e) : = ^ ( a„ + h^Un-A (-l)"e-'"P (6.8a) 

and G2{e) : = ^ [ a„ + ^ T„ff„_ J (-1)"£-'"p, (6.8b) 

where arc the Fourier coefficients of the multiplicatively periodic function h{e) 
discussed just before the statement of Theorem 3.1, and a„, 6„, (T„, and t„ are the 
complex numbers given by 

^_33/2 33/2 35/2 1 

+ — — 7 - „ . . — T + i:K, 



23(Z) + mp) 23(£)-l + mp) 25(D-2 + mp) 

(2m)! (32™+! - 4) 



42'"+i(m!)2(4m2 - l)(32™+i - 2) (D - 2m - 1 + mp) ' 
log 3 (I + 2V3) and (6.9) 



(2m)! (32™+! - 1) 



" ^1 42"-l(™0^(4"^2 - 1)(32™+1 - 2)(-2m - 1 + mp) ' 

In the definition of C7„ we use the Kronecker delta to indicate a term that only 
appears in o-q. 

Remark 6.5. As mentioned in Remark 5.2, the exterior e-neighbourhood of the 

Sicrpinski gasket is the disjoint union of the inner e-ncighbourhood of the tiling and 
the exterior e-neighbourhood of the largest triangle; see Figure 10. This means that 
one can immediately obtain the tube formula for the exterior e-neighbourhood of 
the gasket by adding the tube formula for the tiling SQ and Stcincr's tube formula 
for its convex hull, the equilateral triangle Co (as labeled in Figure 8). However, 
things do not always work out so neatly, as the example of the Koch tiling shows; 
sec Figure 11. In the forthcoming paper [PeWi], precise conditions are given for 
equality to hold as in (6.10). In particular, 

vol<i(r/-*) = Vr(e) + volrf([i^]) (6.10) 

holds whenever one of the following equivalent conditions is verified: 



This set of complex dimensions was predicted in [La-vFl], §10.3, except for the dimensions 
above 0. 
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Figure 10. The exterior e-neighbourhood of the Sicrpinski gasket 
is the disjoint union of the exterior e-neighbourhood of C = [F] 
and the inner e-neighbourhood of the Sierpinski gasket tihng. 




Figure 11. The exterior e-neighbourhood of the Koch curve is 
not simply related to the exterior e-neighbourhood of C = [K] and 
the inner e-neighbourhood of the Koch tiling /C. 




Figure 12. The pentagasket tihng V. 

(1) dT:^F. 

(2) dGq CFioTallq^l,...,Q. 

(3) 5(C\$(C)) C F. 

(4) dC C F. 

(5) _F|5^* DC — for some (and, equivalently, all) e > 0. 

(6) F^^^* = C|"=* \ C for some (and, equivalently, all) e > 0. 

Thus (4.8) (in conjunction with (6.10)) allows one to compute explicit tube formulas 
for a large family of self-similar sets. 

Example 6.6 (The pentagasket tiling). The pentagasket tiling V is constructed 
via the self-similar system defined by the five maps 

= ^a;-|-p„ j-l,...,5, 

with common scaling ratio r = (/)~^, where (/> = (! + •\/5)/2 is the golden ratio, and 
the points = Cj form the vertices of a regular pentagon of side length 1; see 
Fig. 12. 
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The pcntagasket V provides an example of multiple (noncongrucnt) generators 
Gq with q = 1, 2, . . . , 6. Specifically, Gi is a regular pentagon and G2, ■ ■ ■ ,Ge are 
congruent isosceles triangles, as seen in Ti of Fig. 12. This example is developed 
fully in [LaPe2, §9.4]. 
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